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ABSTRACT 


Observations  of  Currents  in  the  Arctic  Ocean 
Robert  L.  Bernstein 

An  array  of  current  meters  was  installed  and  operated  for  three 
months  at  Ice  Island  T-3>  at  85°N,  95°W  in  the  Arctic  Ocean.  Current 
meters  were  placed  at  three  sites  at  40m  and  7Cta  depths,  with  signal 
outputs  all  wired  to  a  common  recording  system.  The  drift  motion  of 
T-3  can  be  eliminated  from  the  records  using  frequent  precise  satellite 
navigation  fixes.  All  records  show  predominant  12  hour  period  inertial 
currents,  superimposed  on  lower  frequency  transient  currents,  which 
occur  during  intervals  of  rapid  ice  drift.  The  transient  currents,  at 
^Om  in  an  upper  layer  veer  to  the  right  of  the  ice  drift  direction, 
while  currents  in  the  lower  layer  at  70m  veer  to  the  left.  The  observed 
ice  drift  and  currents  are  consistent  with  the  theory  of  elementary 
currents  first  advanced  by  Ekman,  and  with  a  recent  theory  by  Pollard. 

The  inertial  currents  have  constant  phase  and  amplitude  which  only 
change  during  intervals  of  ice  displacement,  suggesting  the  local 
generation  of  these  currents  in  the  upper  levels  of  the  ocean. 

The  kinetic  energy  power  spectrum  derived  from  the  current  meter 
time  series  shows  an  even  fall-off  of  energy  proportional  to  the  -5/3 
power  of  frequency.  This  fall-off  is  interrupted  by  an  energy  peak  at 
the  12  hour  inertial  period.  The  -5/3  fall-off  suggests  that  energy 
is  fed  in  at  large  scales  and  long  periods  (>2-3  days)  and  dissipates 
in  an  energy  cascade  through  smaller  scales  and  high  frequencies. 

The  magnitude  of  the  surface  stress  which  the  moving  ice  cover 
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exerts  on  the  underlying  ocean  has  been  estimated  by  a  new  method. 
From  the  decay  of  transient  currents  under  a  non-moving  ice  cover, 
a  simple  theory  gives  a  stress  value  around  0.7  dynes/cm  ,  for 
10  cm/sec  of  relative  velocity  between  the  ice  and  the  water  under 


neath . 
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INTRODUCTION 

Introductory  Remarks 

A  complete  understanding  of  ocean  dynamics  requires  detailed 
observations  of  ocean  currents  over  a  wide  range  of  space  and  time 
scales.  Observation  of  transient  currents  requires  refined  current 
sensors  mounted  on  stable  platforms .  Recent  developments  of  the 
Savonious  rotor  current  meter  and  the  digital  magnetic  tape  recorder 
provide  highly  suitable  instrumentation  for  this  purpose.  The  drift¬ 
ing  ice  stations  on  the  Arctic  ocean  provide  stable  platforms  for 
mounting  the  current  meters  on  rigid  inverted  masts.  The  present  study 
describes  measurements  of  transient  ocean  currents  with  a  three- 
dimensional  array  of  current  meters  with  horizontal  aperture  of  370m 
and  vertical  aperture  of  70m.  The  data  were  all  recorded  on  a  single 
magnetic  tape  providing  a  common  time  base  not  available  in  any 
previous  deep  ocean  current  observations.  The  common  time  base  allows 
precise  comparisons  between  instruments.  The  stable  ice  mounting 
avoids  the  spurious  oscillations  and  noise  injected  by  buoy  mountings 
in  other  oceans.  These  data  are  therefore  unique  in  their  phase 
precision  and  lack  of  high-frequency  noise.  Wave  motion  and  cable 
oscillation  have  been  eliminated  in  this  experiment. 

Historical  Background 

Scientific  observations  of  wind  and  weather,  ice  movement,  and 
the  underlying  Arctic  ocean  dates  back  to  the  earliest  explorers  of 
this  harsh  environment .  The  historic  drift  of  the  "Fram"  by  NANSEN 
(1902)  yielded  results  which  stand  in  many  ways  unimproved  upon. 
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Other  similar  expeditions  followed  Nansen. 

Beginning  in  the  1930’s  Soviet  workers  began  long-term  occupations 
of  drifting  ice  floes  which  continue  to  the  present.  Following  the 
end  of  World  War  II,  U.S.  Air  Force  and  Navy  interests  resulted  in  the 
occupation  of  a  number  of  ice  stations.  Most  permanent  of  these  is 
Fletcher's  Ice  Island  (T-3).  Unlike  ice  floes  which  are  generally 
only  a  few  meters  thick,  T-3  is  many  times  thicker.  It  is  a  piece  of 
glacial  ice,  probably  broken  off  from  Ellesmere  Island.  The  extra 
thickness  keeps  the  ice  intact  so  that  permanent  stations  can  be  main¬ 
tained.  T-3  drifts  in  a  large,  geostrophic  clockwise  circulation, 
called  the  Beaufort  Gyre,  contained  in  the  Canadian  Basin  of  the 
Arctic  Ocean  (figure  l).  T-3  has  been  almost  continuously  manned 

since  1950;  and  during  this  time  it  has  completed  nearly  three  circuits 
around  this  gyre. 

The  drift  of  T-3  and  the  rest  of  the  permanent  ice  cover  in  the 
Canadian  Basin  is  related  to  the  mean  annual  atmospheric  pressure 
distribution  of  the  region. 

General  Drift  of  the  Arctic  Pack  Ice 

A  high  pressure  system  (in  the  Northern  Hemisphere)  is  associated 
with  a  geostrophic  wind  system  causing  clockwise  circulation  around 
the  center  of  the  high.  The  direction  toward  which  the  wind  blows  is 
generally  almost  parallel  to  the  pressure  isobars,  except  near  the 
surface.  There,  friction  and  the  Coriolis  force  interact  in  what  is 
termed  an  Ekman  boundary  layer.  The  resulting  wind  near  the  surface 
blows  outward  from  the  high,  at  an  angle  approximately  30°  to  the 
direction  of  the  isobars  (figure  2a).  For  a  low  pressure  system 
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(figure  2b),  the  geostrophic  circulation  is  counter-clockwise,  with 
surface  wind  blowing  in  toward  the  center  of  the  low. 

In  the  Arctic,  ice  usually  drifts' at  some  acute  angle  (30°  to  45°) 
to  the  right  of  the  direction  of  the  surface  wind  vector  Vg  (figure 
3),  as  was  first  observed  by  Nansen.  From  many  observations  of  Soviet 
ice  stations,  ZUBOV  (1945)  noted  the  resulting  parallel  alignment  of 
pressure  isobars  and  ice  drift  tracks.  FELZENBAUM  (1958)  developed  a 
theory  yielding  parallel  alignments  by  considering  a  model  incorporating 
wind  and  water  stress  forces  on  the  ice,  along  with  Coriolis  force 
(independent  of  latitude),  and  pressure  gradient  force,  while  neglect¬ 
ing  any  large  scale  resistance  to  strain  within  the  ice  cover. 

Zubov’s  result  has  been  considered  valid  for  monthly  and  longer 
time  scales.  Recent  advance  in  navigation  has  allowed  improved 
resolution  of  ice  drift  tracks.  Results  from  the  experiment  to  be 
described  here  showed  that  the  drift  track  generally  conformed  to  the 
isobar  pattern  down  to  a  one  day  time  scale. 

From  Felzenbaum’s  fifteen  year  mean  annual  sea  level  pressure  field 
over  the  Western  Arctic  (figure  4),  CAMPBELL  (1965)  has  calculated 
surface  wind  stress.  He  then  calculated  ice  drift  from  this  stress 
field,  incorporating  the  effects  of  air  and  water  stress,  actual 
Coriolis  variation  with  latitude,  internal  ice  resistance  to  strain, 
and  pressure  gradient  forces.  The  resulting  ice  drift  was  not  quite 
parallel  to  the  isobars,  and  would  produce  a  convergence  of  ice  in  the 
center  of  the  Beaufort  Gyre.  One  of  his  typical  cases  is  shown  in 
figure  5.  The  existence  of  suce  a  convergence  has  yet  to  be  demonstra¬ 
ted.  COACHMAN  (1969)  has  illustrated  that  the  drift  tracks  of  numbers 
of  ice  stations  in  the  Beaufort  Sea  do  roughly  parallel  the  isobars  of 
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T-3  DRIFT  TRACK 


Fig.  1.  Drift  track  of  T-3  from  1952  to  1969* 
covered  in  progressive  vector  diagram,  figure  17 


Insert  shows  area 
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Fig.  2.  (a)  Clockwise  geostrophic  flow  around  a  high-pressure 

system,  with  surface  wind  V$  blowing  outward  30°  from  the  isobars, 
(b)  Counter-clockwise  flow  around  a  low-pressure  system,  V$  blowing 
inward  30°  from  the  isobars . 


Fig.  3*  Ice  drifts  approximately  30°  to  the  right  (in  Northern 
Hemisphere)  of  the  surface  wind  ,  resulting  in  near-parallel 
alignment  of  ice  drift  track  and  isobars . 
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Fig.  k.  Fifteen-year  mean  annual  sea  level  pressure  field  over  the 
Western  Arctic.  From  CAMPBELL  (1965),  after  FELZENBAUM  (1958). 
Pressure  is  in  millibars.  Insert  shows  area  covered  by  figure  17. 
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Fig.  5*  Theoretical  ice  drift  in  the  Arctic,  incorporating  effects 
of  air  and  water  stress,  actual  Coriolis.. variation  with  latitude, 
internal  ice  resistance  to  strain,  and  pressure  gradient  forces. 

From  CAMPBELL  (19 65).  Solid  lines  are  ice  velocity  streamlines;  dash 
lines  are  isotachs  in  cm/sec.  Insert  shows  area  covered  by  figure  17. 
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Fig.  6.  Drift  tracks  of  several  ice  stations  in  the  Beaufort  Sea 
trace  out  the  isobars  of  the  mean  anticyclonic  system  over  the  region. 
From  COACHMAN  and  BARNES  (1961).  Insert  shows  area  covered  by 
figure  17. 


13 


<*t 

23.5 

24.5 

25.5 

26.5 

27.5 

s 

29 

31 

33 

35 

37% 

T 

-2 

-1 

0 

l 

0 

» 

CJ 

Fig.  7.  Typical  Arctic  Ocean  profile  of  temperature  (  T  ),  salinity 
(  S  ),  and  density  (  0^  ).  Note  the  highly  stratified  layer  above 
250m,  atop  a  nearly  uniform  layer.  The  upper  50m  (upper  mixed  layer) 
has  uniform  density  resulting  from  surface  action.  (K.  HUNKINS, 
personal  communication) . 
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Fig.  8.  Dynamic  topography  of  the  Western  Arctic  Ocean,  showing 
convergence  of  surface  water  in  the  center  of  the  Beaufort  Gyre. 
From  COACHMAN  and  BARNES  (1961).  Insert  shows  area  covered  by 
figure  17. 
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the  mean  anticyclonic  system  centered  at  80°  N,  lkO°  W  (figure  6). 

We  now  turn  attention  from  the  ice  cover  to  the  ocean  underneath. 
The  top  50m  of  the  Arctic  Ocean  is  uniform,  being  well-mixed  by 
surface  action,  and  is  referred  to  as  the  upper  mixed  layer  (figure  7). 
Below  this,  a  continuously  stratified  layer  extends  to  a  depth  of  250m, 
which  in  turn  overlies  a  nearly  uniform  lower  layer  extending  to  the 
bottom.  COACHMAN  and  BARNES  (1961)  have  examined  hydrographic  infor¬ 
mation  collected  from  many  sources  over  many  years  from  the  Western 
Arctic.  The  dynamic  topography  which  they  constructed  from  this  data 
(figure  8)  revealed  a  lens  of  surface  water  in  the  center  of  the  gyre, 
associated  with  a  clockwise  geostrophic  gyre  in  these  surface  waters. 
This  gyre  is  quite  similar  to  the  geostrophic  circulation  accompanying 
the  mean  atmospheric  pressure  distribution. 

The  fact  that  the  upper  ocean  geostrophically  follows  the  atmos¬ 
pheric  circulation  may  not  be  too  surprising.  Such  a  result  is  pre¬ 
dictable  from  the  EKMAN  (1905)  theory  of  elementary  currents.  Accord¬ 
ing  to  Ekman,  a  stress  at  the  surface  causes  a  net  transport  of  surface 
waters  at  a  direction  90°  to  the  right  of  the  stress.  For  the  observed 
clockwise  ice  drift  of  the  Beaufort  Gyre,  this  would  result  in  a 
convergence  of  surface  waters  at  the  center  of  the  gyre,  as  seen  by 
Coachman  and  Barnes.  This  convergence  must  be  associated  with  a  clock¬ 
wise  geostrophic  circulation  of  the  water  in  the  Beaufort  Gyre. 

One  of  the  main  results  of  the  present  experiment  is  a  detailed 
description  of  the  manner  in  which  water  converges  toward  the  center  of 
the  gyre.  The  results  will  show  that  transient  currents  are  set  up  by 
storms  moving  the  ice,  which  in  turn  stresses  the  water.  These 
transient  currents,  which  only  last  for  the  duration  of  the  storm,  are 
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shown  to  follow  the  Ekman  theory,  even  on  a  time  scale  as  short  as 
one  day.  Individual  storms  tend  to  move  the  ice  in  apparently  random 
directions,  so  that  it  is  only  over  a  period  of  years  that  the  weak 
Polar  high  pressure  system  manifests  its  presence  with  the  ice  and  also 
the  surface  waters,  causing  net  circulation  around  the  gyre. 
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THE  EXPERIMENT 


The  Current  Meter  Array 

Ice  Island  T-3  provides  a  highly  stable  platform  for  conducting 
oceanographic  studies.  The  island  is  roughly  10km  by  6km  and  30m 
thick,  entirely  surrounded  by  pack  ice  3m  thick  (figure  9).  Other 
advantages  accrue  from  various  independent  operations  on  T-3.  The 
U.S.  Weather  Bureau  collects  wind  and  atmospheric  data,  while  the 
University  of  Washington  has  a  program  involving  hydrographic  casting. 
Lamont -Doherty  Geological  Observatory  supplies  satellite  and  celestial 
navigation,  bottom  and  subbottom  profiling,  gravity  and  magnetics 
measurements . 

In  this  setting,  current  meter  studies  have  been  undertaken  for 
several  years.  The  present  experiment  is  an  attempt  to  improve  or 
extend  previous  results.  The  current  meters  used  were  externally 
recording  (Marine  Advisors,  Model  Q-9)  with  all  telemetry  cabling 
brought  to  a  single  recording  site.  The  Savonius  rotor  speed  sensor 
and  direction  vane  sensor  are  of  a  design  which  has  become  almost 
standard  in  oceanography  for  the  past  decade.  Here  data  were  recorded 
on  a  digital  magnetic  tape  system  (Hewlett-Packard  2012D)  and  a  multi¬ 
point  chart  recorder. 

All  current  meter  outputs  were  sampled  once  a  minute.  This  is 
considered  to  be  high  enough  to  avoid  aliasing  of  speed  channel 
measurements.  Electronic  circuitry  at  the  recording  site  integrates 
speed  measurements  over  a  30-second  interval.  Direction  measurements 
are  subject  to  some  degree  of  aliasing  because  these  measurements  are 
only  integrated  over  a  5-second  interval.  The  resulting  vane  flutter 
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in  the  measurements  is  usually  confined  to  +  5°  or  less. 

Aside  from  a  few  drogue-type  current  measurements  (HUTTKINS,  1966, 
BROWNE  and  CRARY,  1958)*  all  previous  investigators  in  the  Arctic  have 
used  current  meters  with  internal  compasses  for  referencing  the 
direction  of  currents.  Compasses  are  often  unreliable  in  high  lati¬ 
tudes.  In  the  present  experiment,  directions  were  referred  to  the  ice 
itself.  The  current  meters  were  suspended  on  aluminum  pipe  in  an  in¬ 
verted  mast  technique.  The  pipe  rigidly  controls  the  orientation  of 
the  current  meter.  The  orientation  is  determined  once  and  then  assumed 
to  remain  constant . 

All  oceanographic  work  at  T-3  is  done  on  Colby  Bay,  a  stable 
piece  of  pack  ice  permanently  attached  to  the  island  and  much  thinner 
than  the  island  itself.  Current  meters  were  placed  at  depths  of  40m 
and  70m.  The  choice  of  depths  was  based  on  several  considerations. 

From  hydrocast  data  (figure  7)  the  depth  to  the  base  of  the  upper 
mixed  layer  was  determined  to  be  50m.  Since  the  draft  of  the  ice  island 
is  30m,  40m  was  the  logical  choice  for  measurements  in  this  upper  layer, 
in  order  to  minimize  any  effect  of  the  ice  island  wake  on  the  current. 
Previous  experience  with  profiling  the  upper  200m  with  current  meters 
had  revealed  an  interesting  jet-like  current  at  around  70m  (figure  10). 
Thus  the  choice  of  70m,  below  the  upper  mixed  layer,  was  made.  Deeper 
measurements  were  of  course  desirable  but  difficult  to  obtain  because 
of  the  pipe  suspension  system.  One  current  meter  was  placed  at  100m, 
but  failed  to  operate. 

In  earlier  studies  at  T-3,  current  meters  were  placed  in  vertical 
arrays.  In  the  present  experiment,  some  current  meters  were  placed  at 
the  same  depth  but  horizontally  separated  (figure  ll),  in  order  to 
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FLETCHER’S  ICE  ISLAND  (T-3) 


Fig.  9.  Outline  of  Fletcher’s  Ice  Island  (T-3).  The  island  is  a 
30m  thick  piece  of  glacier  ice,  surrounded  by  3m  thick  pack  ice. 

All  oceanographic  work  is  done  on  Colby  Bay,  a  protected  section  of 
pack  ice  attached  to  the  island.  Runway  near  main  camp  is  1.5  km 
long. 
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Fig.  10.  Vertical  profile  of  relative  current  speed  taken  at  T-3 
in  May  1968.  Note  intensification  of  current  centered  around  70m. 
Current  direction  is  toward  l60°T  except  above  30m.  (K.  HUNKINS, 
personal  communication) . 
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Fig.  11.  Current  meter  array  operated  at  T-3  for  the  present  experi¬ 
ment.  (a)  Plan  view  with  Masts  A,  B,  and  C  located  on  Colby  Bay.  For 
scale,  the  runway  shown  on  the  map  is  1.5  km  long,  (b)  Masts  A  and  C 
each  held  two  current  meters  at  40m  and  70m  depths .  Mast  B  held  one 
current  meter  at  40m  depth. 
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TABLE  1  Current  Meter  Characterics 


Site 

Current  Meter 

Depth  (m) 

Days  of  reliable 
record  (from 

0000Z  25  Nov. 

1970 

Comments 

A 

A70 

70 

87 

A40  . 

40 

70 

B 

B40 

40 

Direction  vane  sticks -- 

only  reliable  in  currents 
above  10  cm  sec 

C 

:  C70 

70 

87 

C70  and  C40  appear 
slightly  mare  reliable 
than  A70,  A40  in  small 
currents 

;  C40 

40 

70 

■ 
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examine  horizontal  versus  vertical  coherence  of  currents.  The  main 
characteristics  of  the  various  current  meters  and  the  records  obtained 
from  them  are  listed  in  Table  1.  Three  separate  masts  were  installed 
at  sites  A,  B,  C ,  with  separations  of  85m,  295m,  and  370m  between 
masts.  The  locations  of  sites  A  and  B  were  dictated  by  operational 
constraints.  Site  C  was  placed  as  far  from  A  and  B  as  possible,  again 
limited  by  practical  considerations.  Unfortunately  the  current  meter 
at  site  B  developed  a  sticky  vane  and  was  not  useful  when  current 
speeds  were  much  less  than  lOcm/sec.  Certain  short  stretches  of  the 
record  from  current  meter  B^-0  are  useful  and  were  examined  for  high 
frequency  structure.  The  main  advantage  however,  in  having  records 
at  sites  A  and  C  from  the  same  depths  (40m,  70m)  was  improved  confi¬ 
dence  in  the  working  and  reliability  of  any  one  current  meter.  Later, 
when  current  meters  began  to  fail,  duplication  served  to  extend  the 
useful  total  length  of  the  record. 

Navigation 

Most  reported  current  observations  in  the  Arctic  were  conducted 
prior  to  the  introduction  of  the  Transit  satellite  navigation  system. 
From  the  moving  platform  of  a  drifting  ice  station,  current  meters  are 
only  measuring  currents  relative  to  the  ice  station.  In  order  to 
convert  a  relative  current  vector  to  an  absolute  current  vector  (rela¬ 
tive  to  a  fixed  point  on  the  ocean  bottom,  for  example),  one  must 
vectorially  add  the  ice  platform’s  drift  velocity  to  the  relative 
current  vector.  The  determination  of  ice  drift  velocity  for  this 
conversion  places  great  demands  on  navigation  fix  accuracy,  and  on  the 
frequency  with  which  fixes  are  made.  Using  only  celestial  navigation 
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it  was  often  impossible  for  HUNKINS  (1967)  to  say  whether  currents 
observed  at  ice  island  T-3  were  the  result  of  absolute  currents,  or 
of  island  motions  towing  current  meters  through  the  water.  The  present 
experiment  has  the  advantage  of  frequent,  accurate  satellite  navigation 
fixes.  Thus  the  conversion  of  the  relative  current  velocities  into 
absolute  velocities  becomes  possible.  Fixes  generally  were  obtained 
every  two  hours,  accurate  to  better  than  400m. 
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RESULTS  OF  DATA  ANALYSIS 


Ice  Drift  and  Wind 

As  stated  previously,  the  ice  cover  drifts  mainly  under  the 
influence  of  the  surface  wind,  following  Nansen's  rule.  The  ratio  of 
ice  drift  speed  to  wind  speed  (as  measured  at  10m  elevation)  is 
approximately  1:50.  The  ice  drifts  in  a  direction  about  30°  to  the 
right  of  the  direction  toward  which  the  wind  blows.  To  show  the  degree 
to  which  T-3's  drift  followed  Nansen's  rule  during  the  present  experi¬ 
ment,  the  hourly  averaged  wind  vectors  have  been  rotated  clockwise 
through  30°  and  divided  in  magnitude  by  a  factor  of  50.  The  east  and 
north  components  of  this  vector  is  compared  with  the  two  components 
of  the  ice  island  drift  velocity  in  figure  12.  The  ice  drift  velocity 
was  derived  from  the  satellite  navigation  fixes,  and  filtered  to  suppress 
the  noise  associated  with  fix  inaccuracies.  Variations  in  ice  velocity 
with  periods  less  than  6  hours  were  filtered  out. 

Representation  of  Currents  by  Components 

Ocean  currents  are  vector  quantities  which  change  with  time. 
Horizontal  currents  relative  to  a  moving  platform  were  recorded  at  T-3» 
The  relative  current  vectors  must  first  be  corrected  for  the  motion 
of  T-3  by  the  vectorial  subtraction  of  T-3's  drift. 

A  horizontal  current  vector  may  be  described  by  speed  and 
direction  or  else  by  components,  north  and  east .  Components  are  more 
convenient  to  work  with  when  vector  addition  and  subtraction  are 
required.  All  current  speed  and  direction  records  were  converted  to 
east  and  north  velocity  components,  with  components  linearly  inter- 
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Fig.  12.  Comparison  of  wind  velocity  (solid  line)  and  ice  drift 
velocity  (dash  line).  Wind  velocity  vector  (pointing  in  direction 
toward  which  wind  blows)  has  been  rotated  30°  clockwise  and  divided 
in  magnitude  by  a  factor  of  50  before  taking  east  and  north  compo¬ 
nents.  Ice  drift  velocity  derived  from  satellite  navigation  fixes. 
Nansen’s  rule  is  followed  where  the  two  curves  coincide. 
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Fig.  13.  Relative  currents  (solid  line)  from  ^-Om  current  meter  on 
Mast  C.  Upper  graph  gives  east  component  of  current,  lower  graph 
gives  north  component.  Negative  of  T-3's  velocity  (dash  line)  also 
shown.  Departure  of  solid  line  from  dash  line  gives  absolute  current. 
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TIME  IN  DAYS  FROM  OOOOZ  E5  NOV  1970 


Fig.  Ik.  Same  as  figure  13,  but  for  70m  current  meter  on  Mast  C 
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polated  for  intervals  when  the  speed  fell  below  the  Savonious  rotor 
threshold  speed  of  1  cm/sec.  Since  the  records  contain  information 
about  all  frequencies  from  zero  to  30  cycles  per  hour  (cph),  it  is 
convenient  to  discuss  them  in  terms  of  two  arbitrary  frequency  ranges: 
a  high  frequency  range  of  one  to  thirty  cph,  and  a  low  frequency 
range  of  zero  to  1  cph.  The  low  frequency  range  will  be  discussed 
first.  The  current  meter  records  were  passed  through  a  low-pass 
filter  to  eliminate  frequencies  above  1  cph.  The  resulting  records 
for  site  C,  40m  and  70m  depths  (c40,  C70)  are  shown  in  figures  13  and 
14  respectively. 

The  records  C40  and  C70  shown  in  figures  13  and  Ik  give  relative 
currents  (relative  to  the  drifting  ice  station) .  Also  plotted  on  the 
same  figure  is  the  negative  of  the  ice  velocity.  The  negative  of  the 
ice  velocity  has  been  plotted  to  show  what  the  current  meters  would 
have  recorded,  if  they  had  been  dragged  through  still  water.  In 
figures  13  and  Ik  the  absolute  current  may  be  seen  as  the  departure  of 
the  relative  current  curve  from  the  ice  velocity  curve. 

The  current  meter  records  at  40m  (figure  13)  are  generally 
characterized  by  large  sinusoidal  oscillations  with  a  constant  period 
of  about  twelve  hours,  plus  more  random  low  frequency  behavior.  The 
sinusoidal  oscillations  are  only  prevalent  near  the  beginning  of  the 
record.  The  current  meter  records  at  70ni  (figure  l4)  are  character¬ 
ized  by  persistent  sinusoidal  oscillations  modulated  by  lower 
frequencies.  In  a  later  section  these  strikingly  uniform  sinusoidal 
oscillations  will  be  shown  to  be  inertial  motions.  The  ubiquity  of 
inertial  motion  in  the  ocean  has  only  very  recently  come  to  be  realized. 

Inertial  current  motion  is  a  solution  to  the  simple  linearized 
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equations  of  fluid  motion  in  a  rotating  system  such  as  the  earth. 

By  including  only  the  mass-acceleration  and  Coriolis  terms,  the  equa¬ 
tions  have  a  form 


which  is  satisfied  by 


U  =  A  si*  ft 
\r  =  Ac»sft 


(0 


(2) 


where  (u.vj  are  the  two  components  of  fluid  velocity,  and  £  is  the 
Coriolis  parameter.  On  the  earth,  depends  on  the  latitude  (p  , 
following 

f  =  2  id  sinip  (3) 


where 


/I  =  magnitude  of  the  earth’s  angular  velocity  vector. 


Note  that  the  (  U,V*  )  solution  is  a  circular  rotating  vector,  rotating 
clockwise  with  a  frequency,  £  ,  usually  referred  to  as  the  inertial 
frequency.  At  T-3's  latitude  of  85°30'  N,  the  inertial  frequency 
corresponds  to  a  period  of  12.001  hours,  which  is  just  what  is  observed 
in  the  oscillations  of  figures  13  and  l4.  The  inertial  motion  charac¬ 
teristics  will  be  examined  later  in  more  detail.  First,  the  transient 
currents  with  frequencies  lower  than  the  inertial  frequency  will  be 
examined . 

There  are  major  differences  between  the  current  behavior  at  40m 
in  the  upper  mixed  layer  (figure  13 )  and  at  70®  in  the  continuously 
stratified  layer  below  it  (figure  14).  With  the  exception  of  the  first 
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ten  days  of  the  record,  all  currents  at  ^Om  damp  out  within  a  day  or 
two  after  the  ice  cover  has  ceased  moving.  The  uniform  upper  layer 
apparently  cannot  sustain  a  velocity  shear.  The  inertial  currents  at 
70m  on  the  other  hand  show  little  sign  of  damping  during  intervals 
between  ice  motion  events. 

The  exception  of  no  damping  occurring  during  the  first  ten  days 
of  data  at  ^-Om  might  be  attributed  to  persistence  through  late  autumn 
of  some  stratification  just  under  the  ice,  a  consequence  of  the 
presence  of  summer  melt  water.  Hydrographic  data  collected  through 
the  year  suggests  this  possibility. 

To  examine  the  lower  frequency  transient  current,  the  current 
records  shown  in  figures  13  and  l4  were  passed  through  a  second  low 
pass  filter  to  eliminate  frequencies  at  or  above  the  inertial  frequency. 
Figures  15  and  1 6  show  the  resulting  mean  currents,  here  converted 
from  relative  to  absolute  currents.  The  fixes  from  satellite  navi¬ 
gation  are  placed  on  the  same  figure.  There  is  a  marked  tendency  for 
the  ice  motion  to  impart  a  steady  current  to  the  water,  in  the  same 
direction  as  the  ice  movement.  This  tendency  is  not  as  clear  for  ^-Om 
as  it  is  for  70m,  because  of  the  rapid  damping  in  the  upper  homo¬ 
geneous  layer . 

Examination  of  figure  1 6  shows  some  positive  correlation  between 
ice  displacement  dice  and  the  mean  current  at  7°m  depth,  but  for 
the  north  component  of  motion  only.  The  peculiarly  different  nature 
of  the  first  ten  days  of  the  record,  due  to  persistence  of  summer 
stratification  of  the  upper  layer,  also  shows  up  in  figure  1 6.  The 
mean  current  and  ice  displacement  show  marked  negative  correlation 


during  these  ten  days. 


. 


32 


E 


N 


TIME  IN  DAYS  FROM  OOOOZ  25  NOV  1970 


Fig.  15.  Mast  C,  4Qm  depth  current  meter  records  of  figure  13 
filtered  to  remove  frequencies  above  the  inertial  frequency,  and 
converted  from  relative  to  absolute  current  (solid  line,  scale  to 
left).  Upper  and  lower  graphs  give  East  and  North  components, 
respectively.  Longitude  and  latitude  of  satellite  fixes  (crosses) 
give  components  of  ice  displacement  (scale  to  right). 
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Fig.  1 6.  As  in  figure  15,  but  for  Mast  C,  7Qtn  depth  current  meter. 
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No  clear  correlation  between  and  Vfl  is  evident  in  the  east 

component  of  motion.  Especially  lacking  is  an  east  component  current 
response  to  the  large  ice  displacements  of  days  52  and  59*  The 
reason  for  this  will  become  clearer  in  other  graphic  presentations  of 
the  data  to  be  given  shortly. 

Figure  l6  also  shows  that  any  ice  displacement,  and  any  mean 
current  over  a  period  of  a  month  or  two  was  the  vector  sum  of  transient 
currents  induced  by  a  few  discrete  displacements  of  the  ice  cover. 

These  transient  currents  can  move  in  any  direction,  and  the  addition 
of  one  or  two  ice  displacement  events  to  any  given  month  would  have 
changed  the  monthly  mean  current  significantly. 

Other  Methods  for  Graphical  Presentation  of  Currents 

There  are  always  difficulties  in  graphic  presentation  of  vectors 
measured  at  different  points  in  space  and  time.  There  are  various 
methods,  and  each  has  its  own  merits.  The  component  plot  has  already 
been  presented.  It  tends  to  emphasize  the  high  frequencies  over  the 
low  frequencies.  The  other  methods  to  be  used  are: 

1.  The  progressive  vector  diagram 

2.  The  vector  track  plot 

3.  Spectrum  analysis  results 

Each  of  these  three  graphical  methods,  and  the  features  it  best  empha¬ 
sizes  will  now  be  discussed  in  turn. 

The  Progressive  Vector  Diagram 

In  the  progressive  vector  diagram,  all  the  velocity  vectors  of  a 
time  series  are  placed  head  to  tail,  and  the  resulting  displacement  of 
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a  particle  carried  by  these  currents  is  drawn  out.  Figure  17  gives 
the  progressive  vector  displacement  versus  time  for  currents  at  40m 
and  70m,  and  the  drift  track  of  T-3.  The  progressive  vector  diagram 
is  made  from  absolute,  not  relative  currents.  Note  that  most  of  the 
displacements  are  associated  with  a  very  few  major  storms.  The  mean 
velocity  associated  with  all  three  displacement  vectors  over  the  70 
day  record  is  around  2cm/sec.  The  progressive  vector  diagram  also 
shows,  that  to  a  very  rough  approximation  the  ice  island  and  the  mean 
currents,  both  at  ^Qm  and  70m,  tend  to  travel  together. 

The  area  enclosed  by  figure  17  has  been  added  as  a  detail  to  all 
other  figures  showing  the  whole  Arctic  region. 

In  figure  17  the  surface  wind  vectors  have  been  added  to  the  ice 
island  drift  track.  The  vectors  point  in  the  direction  toward  which 
the  wind  blows.  The  ice  nearly  always  moves  to  the  right  of  the  wind 
vector  by  some  acute  angle  which  is  distributed  around  a  peak  at  30°, 
following  Nansen’s  rule.  This  figure  shows  how  high  winds  during 
storms  correlate  with  rapid  drift  of  the  ice  cover.  All  ice  displace, 
ments  shown  in  figures  15  and  1 6  were  caused  by  storms. 

One  major  trouble  with  progressive  vector  diagrams  of  currents 
measured  with  a  current  meter  is  that  the  point  of  observation  is 
constantly  cl  ging  relative  to  a  particular  water  parcel.  The  pro¬ 
gressive  vector  result  is  not  the  same  as  a  Lagrangian  tag  on  a 
particular  parcel  of  water.  With  these  considerations,  figure  17  may 
be  giving  as  good  agreement  as  could  be  expected  between  a  Lagrangian 
ice  drift  track  and  a  semi-Eulerian  derived  current  track. 
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TIME  IN  DAYS  FROM  OOOOZ  25  NOV  1970 


Fig.  17.  Progressive  vector  diagram,  giving  40m  and  70ni  progressive 
vectors.  T-3  drift  track  also  shown,  with  wind  vectors  along  the 
track,  pointing  in  direction  toward  which  wind  blows.  The  area 
covered  in  this  figure  is  given  as  an  insert  in  other  figures  showing 
the  whole  Arctic  region. 
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Fig.  18 .  Vector  track  diagram,  giving  segments  of  T-3's  drift  track. 
40m  and  70m  low  frequency  transient  absolute  current  vectors  are  placed 
along  the  track.  Time  marked  in  days  from  0000Z  25  Nov  1970  along 
drift  tracks.  Outline  of  T-3  oriented  as  during  experiment  (North  is 
to  top  of  diagram),  with  Colby  Bay  (hatched  lines).  70m  flow  always 
is  to  left  of  ice  drift  direction,  while  40m  flows  either  right  or  left . 
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Vector-Track  Diagram 

To  examine  better  the  low  frequency  transient  current  behavior, 
which  is  most  prominent  during  periods  of  ice  motion,  there  is  another 
method  of  graphical  presentation,  which  here  will  be  called  the 
vector-track  plot.  The  daily  mean  of  the  absolute  current  vector  is 

v 

plotted  with  its  origin  at  the  appropriate  position  along  the  ice 
drift  track.  Both  40m  and  70m  mean  current  vectors  are  drawn  along 
the  drift  track. 

In  figure  l8,  the  vector  track  is  shown  for  the  major  periods  of 
rapid  ice  motion.  The  outline  of  T-3  is  shown  in  the  figure,  properly 
oriented  with  respect  to  true  north.  The  current  vectors  generally 
point  forward  in  the  direction  of  ice  motion  and  to  one  side  or 
another  of  the  ice  drift  track.  The  70m  vector  in  almost  every  case 
points  to  the  left  of  the  drift  track.  The  kOm  vector  points  forward, 
but  either  to  the  right  or  left  of  the  drift  track. 

Further  discussion  of  the  low  frequency  transient  currents  will 
be  reserved  for  a  later  chapter. 

Spectrum  and  Cross  Spectrum  Analysis 

One  last  method  for  presenting  the  data  is  by  obtaining  spectrum 
and  cross  spectrum  results.  The  method  of  BLACKMAN  and  TUKEY  (1958) 
is  used.  Because  of  the  extremely  large  number  of  sampled  points  in 
any  one  of  the  current  records  (l/minute  X  70  days  ^  100,000  pts) 

it  was  desirable  to  obtain  the  spectrum  in  two  parts:  high  frequency 
and  low  frequency. 

The  high  frequency  spectrum  was  obtained  by  choosing  an  interval 
of  data  about  one  day  in  length,  at  a  time  where  little  or  no  ice 
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motion  was  occurring  (during  29-30  November  1970) .  The  low  frequency- 
part  was  obtained  by  using  the  entire  70  day  record  of  the  experiment . 

The  data  had  to  be  "pre-whitened"  by  passing  it  through  a  differ¬ 
encing  filter  (JENKINS,  1968)  before  taking  a  spectrum.  This  cuts 
down  the  low  frequency  power  of  the  inertial  and  lower  frequencies 
which  otherwise  would  completely  contaminate  the  spectrum.  The 
resulting  spectrum  for  C40N  is  shown  in  figure  19  after  the  gain  of 
the  difference  filter  has  been  divided  out  of  the  raw  spectrum. 

The  spectrum  for  C40N  in  figure  19,  as  well  as  all  subsequent 
spectral  estimations,  were  performed  with  20  degrees  of  freedom. 

Tables  2a,  b,  c,  from  GRANGER  (1964)  are  used  to  find  the  95$  confi¬ 
dence  limits. 

To  handle  the  low  frequency  spectral  analysis,  the  low-pass 
filtered  records,  which  were  subsampled  every  l/2  hour  were  used 
(figures  13  and  14).  The  records  were  converted  to  absolute  currents 
before  calculating  spectrums .  Care  must  be  taken  to  divide  out  the 
gain  of  the  low  pass-  filter,  and  to  similarly  correct  the  spectrum  for 
the  effect  of  this  and  any  other  filters  applied  to  the  raw  data,  before 
presenting  a  final,  true  spectrum.  The  l/2  hour  subsampled  data  was 
difference  filtered  for  pre-whitening  purposes  for  this  low  frequency 
analysis  also. 

The  spectrum  of  record  C40N  is  given  in  figure  20,  along  with  the 
high  frequency  tail  of  figure  19.  The  inertial  peak  stands  out  strong¬ 
ly  in  the  spectrum  against  a  background  of  energy  falling  off  with  a 
fairly  uniform  slope  of  around  -5/3  or  -2.  There  is  also  an  increase 
of  energy  between  5  and  10  cph,  including  a  peak  at  10  cph  which  is 
not  as  significant  as  the  inertial  peak.  This  high  frequency  increase 
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Degrees 

of 

Freedom 

TABLE  2a  Power  Spectrum  Confidence  Bands 

Exceeded  by  P%  of  all  Values 

P  95%  90%  10% 

5% 

2 

.  05 

.10 

2.  30 

2.  99 

3 

.  12 

.  20 

2.  08 

2.  60 

4 

.18 

.  26 

1.  94 

2.  37 

5 

.  23 

.  32 

1.  85 

2.  21 

6 

.  27 

.  37 

1.  77 

2.  10 

8 

.  34 

.44 

1.  68 

1.  94 
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.  39 
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1.  83 

12 

.43 

.  53 

1.  55 
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15 

.48 
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1.  48 

1.66 

20 

.  54 

.  62 

1.  42 

1.  55 

30 

.  62 

.69 

1.  34 

1.46 

50 

.69 

.  75 

1.  26 

1.  34 

00 

.  77 

.  82 

1.  18 

1.  22 

(From  GRANGER,  1964) 
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TABLE  2b  Coherence  Confidence  Levels 


Degrees 


of 

freedom 

8 

12 

16 

20 

24 

32 

40 

median 

.  454 

.  361 

.  307 

.  272 

.  247 

.  212 

.  190 

90% 

.  732 

.  607 

.  529 

.475 

.  435 

.  377 

.  338 

95% 

.  795 

.  672 

.  590 

.  532 

.488 

.  425 

.  382 

TABLE  Zc  Confidence  Bands  for  Phase  Angle,  in  Degrees 

95%iconfidence  bands 
Coherence 


Degrees 

of 

freedom 

.1 

.  2 

.  3 

.4 

.  5 

.  6 

.  7 

.  8 

.  9 

.  10 

8 

56 

45 

37 

31 

27 

22 

18 

14 

4 

0 

12 

50 

39 

31 

26 

22 

18 

15 

11 

3 

0 

16 

46 

35 

27 

23 

19 

16 

12 

10 

3 

0 

2  0 

43 

32 

26 

21 

17 

14 

11 

9 

2 

0 

24 

41 

30 

23 

19 

16 

13 

10 

8 

2 

0 

32 

36 

26 

20 

17 

14 

11 

9 

7 

2 

0 

40 

33 

23 

18 

15 

12 

10 

8 

6 

2 

0 

(From  GRANGER,  1964) 


POWER  IN  CM/SEC* *2/CPH 


k2 


FREQUENCY  IN  CPH 


Fig.  19*  High  frequency  spectrum,,  for  40m  current.  Spectrum 
obtained  from  one-minute  sampled  data  for  29-30  Nov  1970.,  with  20 
degrees  of  freedom.  Arrows  give  9 5$  confidence  limits.  A  -5/3 
slope  is  drawn  on  this  log-log  diagram. 
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Fig.  20.  Low-frequency  spectrum  for  40m  current,  including  the 
high-frequency  tail  of  figure  19.  Spectrum  based  on  entire  70  day 
record  of  absolute  current.  Degrees  of  freedom  and  confidence 
limits  as  in  figure  19.  Even  fall-off  of  energy  following  -5/3 
slope  interrupted  by  energy  peak  at  the  12  hour  inertial  frequency 
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Same  as  figure  20,  but  for  70^  currents 


is  associated  with  the  local  Vaisala  frequency. 

Hie  power  spectrum  of  record  C70N  is  shown  in  figure  21.  The 
spectrum  shows  the  same  even  fall-off  of  energy  as  in  the  40m  record. 

It  differs  from  the  i+Om  record  in  having  a  more  intense  inertial  peak, 
and  in  lacking  any  peak  in  energy  around  10  cph.  Both  40m  and  70m 
spectrums  appear  to  reach  a  noise  level  of  2x10"^  ( cm/sec )2/cph  at 
20  cph. 

There  is  some  discontinuity  in  the  spectrums  where  the  high  and 
low  frequency  parts  join  near  1  cph.  This  is  especially  true  for  the 

40m  spectrum.  The  reason  why  the  40m  high  frequency  energy  seems  a 

little  high  is  straightforward.  During  much  of  the  time  while  the 
ice  cover  was  not  in  motion,  the  current  at  ^0m  was  zero.  In  select¬ 
ing  a  one-day  interval  for  the  high  frequency  spectral  estimation,  a 
day  with  non-zero  current  was  chosen,  in  order  to  obtain  an  interesting 

result .  Thus  the  energy  at  1  cph  for  this  day  is  high  in  comparison 

with  the  average  energy  over  the  "JO  day  duration  of  the  experiment. 

Since  the  70m  current  did  not  go  to  zero  like  the  40m  current  during 
periods  of  no  ice  movement,  the  one  day  selected  for  the  high  frequency 
spectral  analysis  was  typical  of  the  average  energy  level  at  that  depth. 

The  -5/3  or  -2  slope  of  figures  20  and  21  on  a  log-log  plot 
means  that  kinetic  energy  of  the  current  follows  a  power  law  of 
frequency,  with  the  exponent  of  the  power  law  given  by  the  measured 
slope.  Similar  results  have  been  seen  before  in  current  meter  records 
by  other  workers  (WEBSTER,  1968),  but  the  present  results  extend  over 
a  wider  range  of  frequencies. 

The  -5/3  slope  is  highly  suggestive  of  an  energy-dissipating 
cascade  from  larger  to  smaller  scales  of  motion,  following  the  theory 
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of  KOLMOGOROFF  (l94l)  for  homogeneous  turbulence.  The  energy  sources 
for  such  a  cascade  may  be  the  low  frequency  large  scale  synoptic 
weather  patterns,  and  the  inertial  motion.  A  -2  slope  would  also 
nearly  fit  the  observed  spectrums,  and  various  theories  have  been 
advanced  predicting  this  slope.  For  example,  PHILLIPS  (1971)  obtains 
a  -2  law  for  any  current  meter  moving  through  an  undulating  well¬ 
layered  medium,  with  fluctuations  only  a  result  of  the  meter  occasion¬ 
ally  passing  from  the  steady  flow  regime  of  one  layer  into  the  steady 
flow  of  another  layer.  K.  BRYAN  (personal  communication)  has  shown 
that  a  field  of  internal  waves  at  their  limiting  slopes  also  would 
show  a  -2  energy  law.  The  data  seem  to  favor  slightly  the  -5/3  law. 

The  interpretation  of  the  even  fall-off  of  energy  with  frequency 
has  yet  to  be  settled.  The  spectrums  of  figures  20  and  21  are  offered 
as  additional  evidence  for  future  theoretical  considerations.  The 
spectrums  have  two  important  points  already  mentioned.  The  high 
frequencies  are  not  affected  by  buoy  or  surface  wave  action,  and  the 
low  frequency  motion  of  the  ice  and  water  tend  to  go  together,  so  that 
the  current  meter  tends  to  follow  the  same  water  mass  while  under  the 
influence  of  synoptic  weather  patterns .  This  would  not  be  true  for 
measurements  made  at  a  bottom-anchored  current  meter  system. 

There  are  a  number  of  cross-spectral  estimates  of  special  interest. 
In  what  follows  the  squared  definition  of  coherence  is  used, 
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where  o,  and  0^  are  the  autospectrums  of  the  two  series ,  and  Co 
and  Quad  are  the  Co-  and  Quadrature  spectrums  of  the  two  series. 

In  the  various  phase  diagrams  for  pairs  of  time  series,  positive 
phase  means  the  first  series  of  the  pair  leads  the  second,  while 
negative  phase  means  it  lags.  The  possible  error  in  the  phase  is 
directly  related  to  the  associated  level  of  coherence,  and  may  be 
obtained  from  Table  2c . 

The  power  spectrum  of  both  the  east  and  north  components  of  C40 
show  almost  exactly  the  same  power  at  the  inertial  peak.  The  cross¬ 
spectrum  of  the  two  components  (figure  22)  shows  little  significant 
coherence  except  around  the  inertial  frequency.  There,  coherence 
approaches  0.8  with  associated  phase  of  east  lagging  north  by  97° 

+  9°.  One  would  expect  a  90°  phase  difference  between  components  for 
a  circularly  rotating  vector  moving  clockwise.  The  cross  spectrum 
of  the  two  components  of  record  C70  (figure  23)  gives  a  similar  picture, 
with  coherence  of  0.9  and  phase  lag  of  96°  +  2°,  at  the  inertial  fre¬ 
quency  . 

Other  cross-spectrums  of  interest  include  the  pairs  C70N-C40N, 
A40N-C40N,  A7ON-C7ON.  There  is  no  significant  coherence  between  7°m 
and  40m  (figure  24)  except  at  the  inertial  frequency,  where  it  reaches 
.65,  with  70m  lagging  40m  by  60°  +  13° .  This  is  in  marked  contrast 
with  the  high  coherence  and  zero  phase  lag  for  the  pairs  at  the  same 
depths  (figures  25  and  26).  Coherence  remains  high  for  frequencies 
up  to  and  beyond  the  inertial  frequency,  dropping  below  the  significant 
level  at  .2  cph  (corresponding  to  a  period  of  5  hours). 
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Fig.  22.  Phase  and  coherence  diagrams  for  cross -spectrum  of  East 
and  North  Components  of  40m  current ,  obtained  from  entire  "JO  day 
record  of  absolute  current .  9 5$  confidence  level  for  coherence  is 

indicated. 
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Fig.  23.  Same  as  figure  22,  but  for  70m  currents 
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Fig.  2k.  Cross -spectrum  results  between  north  components  of  70m 

and  40m  currents  measured  at  Mast  C.  Only  significant  coherence 
occurs  at  the  inertial  frequency  -f  . 
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Fig.  25.  Same  as  figure  2k,  but  for  currents  measured  at  YOm  from 
Mast  A  and  at  YOm  from  Mast  C.  Horizontal  separation  between  masts 
is  370a1.  Significant  coherence  up  to  .2  cph  (5  hour  period). 
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Fig.  26.  Same  as  figure  25,  but  for  currents  measured  at  the  70m 
level  at  Masts  A  and  C.  Significant  coherence  up  to  .2  cph. 
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THEORY 


The  Pollard  Theory  of  Geostrophic  and  Inertial  Currents 

POLLARD  (1970)  has  summarized  contributions  to  the  theory  of 
inertial  and  geostrophic  current  theory  by  ROSSBY  (1937)>  STOMMEL 
and  VERONIS  (1956)  and  others.  He  presented  a  model  including  effects 
of  Coriolis  force,  pressure  gradients,  and  density  stratification.  To 
model  the  driving  stress  at  the  surface,  due  in  his  case  to  the  wind, 
and  in  the  Arctic  due  to  ice  movement,  he  considered  that  momentum 
injected  into  the  ocean  by  stress  at  the  surface  is  rapidly  and  evenly 
distributed  through  the  upper  mixed  layer  by  mixing  processes.  For  a 
valid  model  the  time  for  this  mixing  must  be  short  compared  with  the 
inertial  period.  The  surface  stress  is  then  replaced  by  a  body  force 
distributed  equally  through  the  upper  mixed  layer. 

Pollard  begins  with  the  linearized  fluid  equations  of  motion  and 
continuity,  which  take  the  form 


Ik-  w=° 
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b* 
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where  U  and  V  are  velocities  in  the  X  and  y  direct ions ,  w  is 
vertical  velocity,  b  is  the  buoyancy  acceleration  sV? 
where  density  and  Mfc)  is  the  Brunt -Vaisala  frequency 

defined  by  ■  f  is  Coriolis  ps*r or  •  is 

th0  body  forco  oquivalont  stross^  and  Zfe)  is  its  distribution  with 
depth . 

A  solution  to  these  equations  may  be  obtained  by  the  method  of 
normal  modes.  One  solution  by  Pollard  (his  case  II ),  was  for  a 
strongly  stratified  layer  of  thickness  overlying  a  weakly  stratified 
bottom  layer  of  thickness  H^.  This  solution  (Equations  6, 7,8)  can  be 
separated  into  oscillatory  and  non-os dilatory  parts.  The  non¬ 
oscillating  part  represents  the  geostrophic  or  mean  current  contri¬ 
bution.  The  oscillating  part  has  eigen  frequencies  St  which  all  lie 
within  a  few  percent  of  the  inertial  frequency  ^  for  a  wide  range  of 
typical  ocean  parameters. 
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Solutions  were  calculated  from  the  Equations  by  digital  computer, 
and  checked  by  reproducing  several  of  Pollard’s  figures.  The 
solution’s  general  characteristics  are  sketched  as  follows  for  a 
stress 


Tc  ~  o 


all  t 

r.fc©*ta<r)f  ^ofoVr 

0  tfftterwiJC. 


In  a  hodograph  plot  of  the  two  velocity  components,  the  non¬ 
oscillating  part  of  the  solution  steadily  increases  while  the  stress 
is  applied  (figure  27a).  During  this  time  the  oscillating  part  of  the 
solution  (  )  follows  the  circle  A  for  currents  in  the  upper 

mixed  layer,  and  circle  B  for  currents  below  the  upper  mixed  layer,  as 
shown  in  figure  27a.  The  oscillatory  part  of  the  solution  reaches  a 
maximum  velocity  for  times  , where  TJ-  is  the  inertial 

period,  and  goes  to  zero  for  f...  .  Note  that  while  the  stress 

is  being  applied,  the  oscillating  part  has  a  mean  transverse  flow 
associated  with  it.  This  flow  is  to  the  right  of  the  stress  direction 
in  the  upper  mixed  layer  and  to  the  left  below  it .  This  transverse 
flow  tilts  the  isopycnal  surfaces  down  to  the  right  of  the  stress 
direction.  This  increasing  tilt  must  then  be  associated  with  the 
steadily  increasing  (non-oscillatory)  geostrophic  current  ^  , 

flowing  in  the  direction  of  stress. 
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When  the  stress  returns  to  zero,  the  non-oscillating  part  of  the 
solution  thereafter  retains  the  value  it  had  attained  at  the  time  T  . 

The  oscillating  part  of  the  solution  changes  to  follow  a  circle  C 
centered  on  the  origin,  with  radius  r(M)  -  (Ux(t>r)  +  V*  ft*T  > ) 1 
(figure  27c).  Hie  time  to  complete  one  orbit  or  circle  C  is  close  to 
but  slightly  less  than  one  inertial  period.  Thus  the  mean  transverse 
flow  has  disappeared  along  with  the  stress,  but  the  isopycnals  remain 
at  their  final  tilt  position,  in  balance  geostrophically  with  a  constant 
current  .  For  times  t>?  ,  the  radius  r  of  circle  C  changes 

slowly  compared  to  the  inertial  period,  and  this  change  is  dependent  on 
the  horizontal  scale  features  of  the  previously  applied  stress  field. 

From  inspection  of  figure  2Jb  for  times  when  stress  is  being 
applied,  current  vectors  averaged  over  one  or  two  inertial  periods 
make  an  acute  angle  with  the  stress  direction.  For  the  solution  in  the 
upper  mixed  layer,  the  averaged  current  vector  veers  to  the  right  of 
the  stress  direction,  while  below  the  upper  mixed  layer,  the  averaged 
current  veers  to  the  left  of  the  stress  direction.  Qualitatively,  this 
character  of  the  solution  may  be  compared  with  the  results  of  the 
present  experiment,  exhibited  in  the  vector-track  diagram  (figure  18). 

Thus  far,  mainly  the  time  dependency  of  the  solution  has  been 
discussed.  The  vertical  profile  of  the  current  depends  on  the  strati¬ 
fication.  For  the  two-stratified-layer  model  (Pollard’s  Case  II ) 

Pollard  has  provided  vertical  profiles  for  typical  deep  ocean  para¬ 
meters.  The  solutions  give  mostly  barotropic  geostrophic  or  mean 
current  through  the  deeper  layer,  with  the  possibility  of  intensifi¬ 
cation  upward  through  the  upper  layer.  The  inertial  current  part  of 
the  solution  also  has  maximum  amplitude  in  the  upper  mixed  layer  and 
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stress  TJ 


Fig.  27.  Hodograph  sketch  of  Pollard  solution  characteristics. 

(a)  Stress  turned  on  ‘f’O  to  *(=7*  }  geostrophic  current  ^  builds 

in  the  direction  of  stress,  while  inertially  oscillating  current 
(  Ut  ,  Vt  )  follows  circles  A  (upper  mixed  layer)  and  B  (below  upper 
mixed  layer),  (b)  Total  current  equals  sum  of  /g  and  (  Ur,  Vj  ), 
for  upper  mixed  layer  (solid  line)  and  below  upper  mixed  layer  (dash 
line),  (c)  After  stress  turns  off,  geostrophic  current  remains 
constant,  while  inertial  current  switches  from  circle  A  to  circle  C. 
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just  below  it,  decreasing  rapidly  with  increasing  depth.  The  inertial 
currents  are  very  small  throughout  the  lower  layer. 

Application  of  the  Pollard  Model  to  the  Arctic 

For  the  next  step,  parameters  more  typical  of  the  Arctic  were 
substituted.  Water  depth  was  taken  as  1500m,  upper  and  lower  layer 
stratifications  chosen  corresponding  to  Vaisala  periods  of  lk  and  70 
minutes,  with  layer  thicknesses  of  250  and  125m  respectively.  An 
upper  mixed  layer  was  taken  as  the  upper  50m  of  the  250m  thick 
stratified  layer.  Surface  stress  is  1  dyne/cm^.  We  are  left  to 
choose  k  ,  where  L,  is  the  horizontal  scale  of  the  surface 

stress  field,  and  also  a  position  X  in  this  field  at  which  to  calculate 
a  solution.  Kx  will  be  taken  such  that  S’rnlkx)*)  for  convenience. 
Figures  28a, b  show  the  resulting  vertical  current  profiles  at  time 
inertial  periods,  where  stress  was  turned  on  from  "t-0  to  5* 
inertial  periods. 

The  vertical  profiles  of  current  from  the  solution  show  different 
results  for  k  corresponding  to  horizontal  ice  stress  scales  of  125  and 
500  km.  There  is  no  a  priori  reason  to  choose  any  particular  k,  but 
the  smaller  scale  gives  better  agreement  with  profiling  observations 
(figure  10). 

Pollard's  model  for  the  wind  stress  over  an  open  ocean  requires 
an  important  modification  for  an  ice-covered  ocean.  In  the  open  ocean, 
when  a  storm  has  passed  and  excited  currents  in  the  upper  mixed  layer, 
these  currents  are  free  to  continue.  In  an  ice-covered  ocean  the  ice 
usually  stops  moving  before  or  shortly  after  passage  of  a  storm. 
Currents  created  in  a  uniform  upper  mixed  layer  cannot  sustain  a 
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Fig.  28.  Vertical  profiles  of  currents  for  Pollard  solution  with 
typical  Arctic  parameters,  (a)  Horizontal  scale  of  stress  L  =  500  km. 
(b)  L  =  125  km. 
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velocity  shear  with  respect  to  this  non-moving  ice  cover,  and  rapidly 
lose  their  energy.  If  this  is  the  case,  the  current  profiles  of 
figure  28b  would  be  slightly  modified,  with  zero  current  above  50m, 
which  is  somewhat  more  in  agreement  with  vertical  current  profile 
measurements,  as  in  figure  10. 

Validity  of  the  Pollard  Model 

One  basic  physical  assumption  of  the  Pollard  theory  requires 
rapid  mixing  of  stress  through  the  upper  mixed  layer.  The  layer,  be¬ 
cause  of  its  uniform  density  distribution,  is  not  supposed  to  be  able 
to  sustain  a  mean  velocity  shear  in  the  vertical.  With  the  assumption 
of  mixing  on  a  time  scale  short  with  respect  to  the  inertial  period, 
the  surface  stress  can  then  be  replaced  by  a  body  force  distributed 
through  the  entire  upper  mixed  layer. 

Evidence  for  rapid  mixing  comes  from  a  Ross  100khz  echo  sounder 
record  operated  near  the  current  meter  array.  The  echo  sounder  often 
showed  a  strong  return  from  planktonic  organisms  floating  at  the  base 
of  the  upper  mixed  layer  (HUNKINS,  197l).>  at  a  depth  of  50ni  when  these 
current  measurements  were  made.  The  depth  to  the  base  of  the  upper 
mixed  layer  thus  can  be  monitored  continuously,  with  an  accuracy  of 
about  0.5m.  Small  undulations  often  occur,  with  amplitudes  of  about 
1m,  periods  of  a  few  minutes  and  longer.  There  are  certain  drawbacks 
in  this  sounder  data.  The  organisms  sometimes  are  not  present  so 
there  is  no  return.  At  other  times  there  are  two  or  more  returns 
separated  by  a  few  meters  and  it  becomes  difficult  to  follow  any  one 
return  for  more  than  a  day. 

Figures  29a,  b,  show  the  Ross  record  at  2  and  k  hours  following 
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onset  of  an  ice  motion  event.  Figure  29a  shows  the  interface  fairly 
quiescent  and  still  typical  of  times  of  no  ice  movement.  Figure  29b, 
k  hours  after  onset  of  ice  movement ,  shows  the  beginning  of  vigorous 
motion  at  and  above  the  ’pOm  interface,  with  little  disturbance  below 
the  interface.  With  this  supporting  evidence  for  rapid  mixing,  the 
Pollard  theory  will  now  be  used  as  a  framework  for  further  discussion. 
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Fig.  29.  Photographs  of  two  one-hour  segments  of  the  Ross  100  khz 
echo  sounder,  2  hours  (upper  photo)  and  k  hours  (lower  photo)  after 
onset  of  ice  movement  at  1100Z  5  Dec  197° •  Base  of  thick  upper 
mixed  layer  clearly  shows,  with  disturbances  occurring  in  lower  photo. 
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DISCUSSION  OF  TRANSIENT  AND  INERTIAL  CURRENTS 


Details  of  Transient  Current  Motion 

It  was  previously  stated  that  the  mean  (or  geostrophic)  absolute 
current  and  the  ice  displacement  have  a  rough  correspondence  as  shown 
in  figures  15  and  16.  The  correlation  of  changes  in  mean  current  with 
ice  displacement  can  be  understood  if 
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where  is  the  change  in  mean  or  geostrophic  current,  is 

«-p 

the  stress  exerted  by  the  ice  on  the  water,  and  I  is  duration  of 
stress.  Equation  10  is  taken  from  POLLARD  (1970),  and  is  just  a 
statement  of  Newton's  second  law  of  motion.  It  is  convenient  to 
assume  that  is  approximately  proportional  to  ice  velocity,  V. 
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although  a  V  law  would  probably  be  more  physically  accurate. 
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where  ditC  is  the  ice  displacement.’ 

The  vector  track  plot  of  the  mean  currents  (figure  l8)  can  be 
interpreted  as  showing  the  effect  of  Ekman  transport  to  the  right  of 
the  stress  direction  above  the  ^0m  interface,  and  to  the  left  below 
the  interface.  The  Ekman  transport  causes  the  isopycnals  to  tilt 
down  90°  to  the  right  of  the  ice  displacement  direction.  This  trans¬ 
port  must  be  combined  with  the  geostrophic  current  which  builds  up  in 
the  direction  of  ice  stress,  as  the  isopycnals  begin  to  slope.  Such 
a  condition  is  shown  in  figure  27b.  The  result  of  these  cross  and 
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down  stream  transports  would  be  a  vector  making  some  acute  angle  with 
the  drift  track.  This  angle  would  be  to  the  right  of  the  ice  displace¬ 
ment  direction  for  currents  within  the  upper  mixed  layer ,  and  to  the 
left  below  this  layer. 

Ekman  type  transport  has  previously  been  observed  from  an  ice  floe 
by  HUWKINS  (1966).  The  ice  floe  had  the  same  thickness  as  the  surround¬ 
ing  pack  ice,  unlike  the  anomalous  30m  draught  of  T-3.  An  Ekman  spiral 
causing  net  transport  to  the  right  of  the  direction  of  ice  movement 
was  measured,  with  most  of  the  transport  confined  to  the  top  l8m  of 
the  upper  mixed  layer.  Mean  flow  to  the  left  at  70m  would  be  the 
compensating  current  to  satisfy,  conservation  of  mass  under  these 
conditions . 

In  the  present  experiment,  during  about  half  of  the  cases  in  the 
vector-track  plots  of  figure  l8,  the  ^Qm  vector  joins  the  7Qdi  vector 
on  the  left  side  of  the  drift  track.  One  possible  explanation  for  this 
could  be  that  these  near  surface  measurements  are  being  made  too  close 
to  the  draught  of .T-3,  which  disturbs  them.  This  possibility  may  be 
checked  by  examining  the  4-Gm  mean  flow  on  the  vector  track  plot, 
along  with  the  orientation  of  T-3  during  the  experiment.  When  the 
ice  island  moves  in  such  a  way  as  to  keep  the  current  meters  ahead  of 
the  island,  one  might  expect  the  minimum  effect  of  the  island's  draught. 
Inspection  of  figure  18  shows  the  current  meters  more  at  the  "bow" 
of  T-3  for  intervals  17-19,  19-22,  31-3^  38-^0,  and  more  at  the 
"stern"  for  intervals  10-12,  13-15,  ^3-^5,  59-81,  66-67.  The  case 
51-55  is  marginally  in  between.  Note  that  the  first  four  cases  have 
developed  kOm.  flow  to  the  right  of  the  drift  track,  while  the  latter 


five  have  not . 
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A  second  possible  explanation  for  the  uncertain  directional 
behavior  of  the  40m  flow  could  be  from  the  change  in  depth  of  the  50m 
interface  following  sloping  of  the  isopycnals.  The  Ross  echo  sounder 
record  showed  that  the  interface  came  up  to  about  3Cta  during  the  ice 
movement  (10-12),  and  to  about  4lm  for  the  movement  (13-15)*  The 
vector  track  for  these  two  cases  shows  the  h-Om.  current  to  the  left  in 
the  former  case,  and  indecisive  in  the  latter  case.  A  sufficiently 
shallow  interface  depth  would  put  both  current  meters  in  the  lower 
layer.  The  70m  instrument  was  twice  as  far  from  the  interface  as  the 
40m  instrument.  Therefore  it  is  much  less  likely  that  isopycnal  tilt¬ 
ing  could  deepen  the  interface  to  70m  depth,  and  thus  cause  the  ^Om 
vector  to  point  to  the  right  of  the  drift  track.  This  in  fact  was 
never  observed  to  happen.  The  early  termination  of  the  Ross  echo 
sounder  record  did  not  allow  comparisons  for  later  ice  displacements. 

The  two  abovementioned  possible  explanations  suggested  as  giving 
the  uncertain  directional  behavior  of  the  40m  current  vector  also  may 
explain  why  the  b-Om  current  magnitude  is  less  than  the  70m  current 
magnitude . 

From  the  observations  shown  in  figures  15  and  1 6,  it  is  actually 
more  correct  to  say  that  the  mean  flow  builds  while  the  ice  is  moving, 
but  when  the  ice  halts,  the  mean  flow  rapidly  drops  back  to  zero. 

This  decline  in  mean  current  occurs  both  at  40m  and  70m*  Since  this 
decline  in  mean  flow  occurs  at  70m,  even  though  inertial  currents 
constantly  persist  there,  frictional  dissipation  cannot  be  invoked  as 
the  agent  causing  the  disappearance  of  the  mean  (geostrophic )  current. 
A  frictional  mechanism  would  erode  both  geostrophic  and  inertial 
motion.  No  acceptable  explanation  can  presently  be  offered  for  under- 
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standing  the  flow  situation  following  cessation  of  ice  cover  motion. 

The  full  internal  wave  equations  may  he  involved,  including  effects 
of  density  stratification,  Coriolis  force,  pressure  gradients  and 
viscous  forces . 

Synoptic  Considerations  of  Transient  Currents 

TheOOOOZ  daily  surface  pressure  charts  for  the  Northern  Hemi¬ 
sphere  produced  by  the  U.S.  Weather  Bureau  cover  the  area  of  the  present 
experiment.  For  nearly  every  one  of  the  70  days  of  the  experiment,  the 
ice  island  drifted  parallel  to  the  pressure  isobars,  keeping  higher 
pressures  to  the  right  of  the  drift  track.  This  is  in  accord  with  the 
early  work  of  ZUBOV  ( 19^-5  )>  referred  to  in  the  first  chapter,  and 
known  as  Zubov's  rule. 

The  degree  to  which  the  ice  island  conforms  to  Zubov's  rule  is 
illustrated  for  the  drift  of  days  38  to  ^0,  best  seen  in  the  vector- 
track  figure  l8.  The  pressure  charts  for  the  three  days  are  reproduced 
in  figure  30*  T-3’s  drift  velocity  is  placed  on  the  charts.  There 

was  negligible  drift  at  the  beginning  of  day  38 ,  and  the  pressure  chart 
shows  very  weak  pressure  gradients  at  T-3's  position  then.  On  days 
39  and  40  a  weak  low  moved  by  T-3  going  toward  the  pole,  and  the  island 
drift  velocity  direction  swung  around  to  remain  parallel  to  the 
changing  isobar  pattern. 

In  the  example  of  these  three  days,  the  70m  mean  absolute  current 
always  provides  some  component  of  flow  toward  the  center  of  the  low. 

The  40m  current  had  a  component  outward  from  the  low.  Such  a  result 
is  in  accord  with  the  simple  Ekman  current  system  theory. 

It  is  interesting  to  note  that  a  similar  result  has  been  obtained 
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Fig.  30*  U.S.  Weather  Bureau  daily  sea  level  pressure  charts  for 
days  38^  39*  ^0  (after  25  Nov  1970)  of  the  experiment.  Isobars  are 
drawn  with  ^  millibar  increments.  North  Pole,  90°  W  longitude,  and 
T-3  position  (85°  N,  95°  W)  are  plotted.  T-3's  daily  drift  vector 
shown,  with  numerals  indicating  drift  rate  in  km/day.  Drift  vector 
tends  to  remain  parallel  to  isobars . 
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10  14 

days 


85.3°N 


85.2°N 


Fig.  31 •  Detail  of  70m  absolute  current,  north  component  (solid  line). 
Time  is  in  days  from  0000Z  25  Nov  1970.  Latitude  from  satellite  fixes 
also  plotted  (crosses).  Note  that  inertial  oscillation  pattern  only 
changes  at  times  of  ice  displacement. 
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once  before.  BROWNE  and  CRARY  (195^)  placed  a  drogue-type  current 
meter  at  150m  below  T-3  in  1952,  and  saw  a  flow  at  this  depth  in  a 
direction  toward  lower  pressure  on  weather  charts. 

Details  of  Inertial  Current  Motion 

The  proper  identification  must  be  made  of  the  12  hour  oscillations 
which  are  most  prominent  in  the  70111  currents  (figure  lk) .  The  only 
justification  so  far  offered  for  identifying  these  oscillations  as 
inertial  currents  has  been  that  their  period  is  very  close  to  the  local 
12  hour  inertial  period.  But  are  these  oscillations  inertial  currents 
or  semi-diurnal  tides?  A  semi-diurnal  tidal  current  having  a  period 
of  12  hours  25  minutes  would  be  nearly  indistinguishable  from  an  in¬ 
ertial  current  at  high  latitudes.  Several  pieces  of  evidence  point  to 
the  inertial  character  of  the  oscillations. 

1.  The  intermittent,  transient  nature  of  the  oscillations,  is 
closely  associated  with  variations  in  the  wind  and  ice  motion  (HUNKINS, 
1967).  The  phase  and  amplitude  of  the  oscillations  change  abruptly 
with  onset  of  ice  movement.  This  is  clearly  illustrated  by  the  detail¬ 
ed  examination  of  absolute  current  before,  during,  and  after  a 
particular  ice  displacement  event  (figure  31) •  Thus  oscillations  show 
no  change  in  phase  or  amplitude  until  the  ice  actually  moves. 

2.  The  different  character  of  the  oscillations  in  and  below  the 
upper  mixed  layer  tends  to  eliminate  the  possibility  of  the  barotropic 
tide,  which  would  have  uniform  current  velocity  at  all  depths. 

3.  Free  internal  waves  of  semi-diurnal  period  are  not  allowable 
above  latitude  of  72°,  because  their  frequency  is  then  less  than  the 
local  inertial  frequency.  Free  internal  waves  are  restricted  to  the 
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frequency  band  between  the  inertial  and  the  Vaisala  frequencies. 

k.  Surface  tides  at  Arctic  Ocean  coasts  are  small.  At  Barrow, 
Alaska,  the  tidal  amplitude  is  less  than  15  cm.  Barrow  has  an  off¬ 
shore  shelf  which  would  tend  to  amplify  the  tide.  Tides  in  the  deep 
Arctic  Ocean  may  only  be  on  the  order  of  a  few  centimeters  in  height, 
with  associated  currents  reaching  perhaps  a  small  fraction  of  a  centi¬ 
meter  per  second. 

5 .  Although  the  time  series  length  of  70  days  collected  in  this 
experiment  does  not  allow  sufficient  resolution,  all  power  spectra 
show  the  inertial  peak  centered  more  on  the  local  inertial  frequency, 
than  on  the  semi-diurnal  tidal  frequency. 

The  five  abovementioned  points  may  be  taken  as  strong  evidence 
for  the  inertial  character  of  the  12  hour  oscillations. 

The  dependence  of  inertial  currents  on  the  ice  displacement  as 
illustrated  in  figure  31>  points  to  the  local  generation  at  the  surface 
of  these  currents.  POLLARD  and  MILLARD  (1970)  came  to  a  similar 
conclusion  from  examination  of  current  and  wind  records  from  the 
Atlantic  Ocean. 

The  inertial  currents  observed  on  record  C70  (figure  lA)  remain 
mostly  constant  in  amplitude,  phase  and  frequency  during  intervals  of 
no  ice  movement.  Comparing  the  oscillations  before  and  after  various 
similar  ice  motion  events,  sometimes  they  are  built  up,  sometimes 
almost  destroyed. 

The  Pollard  solution  for  inertial  currents  (Equations  7>8) 
depends  primarily  on  T  >  the  duration  of  surface  stress,  in  compari¬ 
son  with  the  inertial  period.  If  T  equals  the  inertial  period,  the 
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solution  gives  zero  inertial  current  amplitude.  This  can  be  under¬ 
stood  physically  as  follows:  during  the  first  half  of  the  inertial 
period,  inertial  motion  is  built  up  moving  down  wind,  then  the 
current  begins  to  turn,  and  by  time  T/2 ,  has  begun  to  head  upwind. 

If  the  stress  persists  in  the  same  direction  during  the  second  half 
of  the  inertial  period,  the  stress  will  knock  down  any  inertial 
current  built  up  during  the  first  half. 

Stress  durations  equal  to  an  integral  number  of  inertial  periods 
yield  no  inertial  current  if  no  such  current  existed  prior  to  the 
onset  of  the  stress.  Suppose  an  inertial  current  of  a  particular 
phase  and  amplitude  had  existed  prior  to  onset  of  a  stress  lasting  an 
integral  number  of  inertial  periods.  Then  by  linearity  the  inertial 
current  would  return  to  its  initial  phase  and  amplitude  after  the 
stress  was  turned  off.  Time  changes  in  stress  which  are  short 
compared  with  the  inertial  period  are  most  effective  in  producing 
inertial  oscillations. 

The  observable  changes  in  inertial  current  amplitude  may  be 
attributed  to  differences  in  stress  duration  T  ,  but  there  are  dif¬ 
ficulties  in  determining  T  from  satellite  navigation  fixes,  to  the 
accuracy  necessary.  A  difference  of  three  hours  in  the  determination 
of  T1  gives  very  different  inertial  current  amplitudes.  Also, 
certain  small  amplitude  but  short  time  variations  in  ice  movement 
might  have  a  larger  influence  than  a  large  amplitude,  longer  time 
variation,  yet  the  short  time  variations  can  go  undetected  by  the 
navigation  fixes.  For  these  reasons  it  has  been  impossible  to  find  a 
meaningful  relation  between  changes  in  inertial  motion  and  ice  dis¬ 
placement,  except  to  say  that  changes  only  occur  following  ice  movement. 
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Estimation  of  Surface  Stress  from  Transient  Current  Decay 

The  rapid  decay  in  40m  current  speed  following  ice  movement  as 
observed  in  figure  Ik  can  be  used  to  obtain  surface  stress  values. 
The  simplest  possible  model  of  stress  proportional  to  velocity  gives 


a 


where  is  the  current  at  40m,  A  is  a  proportionality  constant, 

r  is  the  surface  stress,  fl  the  thickness  of  the  upper  mixed  layer, 
and  ^  the  density  of  seawater.  Equation  13  has  an  exponential 
solution 


Values  for  A  are  obtained  by  plotting  these  transient  decays  which 
follow  a  number  of  storms  on  semi-log  velocity  versus  time  coordinates 
(figure  32).  The  resulting  slopes  provide  values  of  A  listed  in 
Table  3.  From  Equation  13  T-y Aav^  ,  and  for  the  average  value 

of  A  from  Table  3  and  a  characteristic  10  cm  S€ t 1  we  find 

^  O.GS  dynes/cm2.  This  value  is  in  rough  agreement  with 
ice  stress  values  obtained  by  other  estimates,  including  velocity 
shear  measurements  made  in  the  upper  few  meters  of  water  near  a  moving 
ice  cover  (K.  HUNKINS,  personal  communication). 
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SPEED,  CM/SEC 


40m  Current  Decays 
•  DEC  5 


INERTIAL  PERIODS  FROM  END 
OF  ICE  MOVEMENT 


Fig.  32.  Transient  decay  of  kOm  low-frequency  transient  current 
speed  following  end  of  ice  displacement.  Exponential  character  of 
decay  shown  by  straight  line  fall-off  on  semi-log  coordinates. 


TABLE  3  Slopes  From  Figure  32  Giving  Frictional  Stress  Constant  A 


Event 

Slope 

6 

tan0 

A(sec 

-i)  =  - 

tan(9  *2.3 

8.  6xl04 

Dec. 

5 

25° 

.466 

1.  23 

10’5 

Dec. 

9 

42° 

.  900 

2.  38 

io-5 

Dec. 

14 

30° 

.  577 

1.  53 

10‘5 

Jan. 

15 

16° 

.  287 

.  76 

10~5 

Jan. 

23 

33° 

.649 

1.  72 

1— ■ 

o 

i 

U1 

Jan. 

30 

12° 

.  213 

.  57 

10‘5 

•8.  19xl0~5 
Average  A  =  6 


=  1.  36xl0-5  sec-1  (=  TT-. - - ^— ) 

1.  7  inertial  periods 
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CONCLUSION 


Energy  input  into  the  Arctic  Ocean  upper  levels  is  by  wind  stress 
on  the  ice  cover,  which  in  turn  stresses  the  water.  The  magnitude  of 
this  stress  has  been  estimated  by  a  new  method  to  be  about  0.7 
dyne s/ cm  .  The  surface  stress  varies  in  time  in  the  same  fashion  as 
the  synoptic  weather  scale.  The  Arctic  Ocean  responds  to  this  stress 
with  transient  currents,  primarily  in  the  baroclinic  mode.  From 
current  meter  observations,  the  upper  mixed  layer  moves  often  to  the 
right,  but  sometimes  to  the  left  of  the  direction  of  ice  displacement. 

A  deeper  layer,  characterized  by  the  70m  level,  moves  universally  to 
the  left  of  the  stress  direction.  These  observations  are  qualita¬ 
tively  in  agreement  with  the  elementary  current  system  advanced  by 
Ekman,  and  with  a  recent  theory  by  Pollard. 

The  kinetic  energy  power  spectrum,  derived  from  current  meter 
time  series  measurements,  indicates  that  energy  is  fed  in  at  some  long 
period,  T  2-3  days.  This  energy  probably  dissipates  in  an  energy 
cascade  through  higher  frequencies,  as  an  even  fall-off  of  energy 
proportional  to  the  -5/3  power  of  frequency.  The  even  fall-off  is 
interrupted  by  an  energy  peak  at  the  12  hour  period,  corresponding  to 
the  local  inertial  period.  The  presence  of  strong  currents  of 
inertial  period  is  a  resonance  phenomenon  associated  with  variable 
ocean  currents  measured  on  the  rotating  earth.  These  inertial  currents 
have  the  form  of  an  unmodulated  sinusoid  which  is  only  significantly 
disturbed  during  the  occurrence  of  local  ice  movement.  Inertial 
currents  in  the  upper  levels  of  the  Arctic  Ocean  are  therefore  locally 
generated.  Attempts  at  more  quantitative  comparison  of  inertial  current 
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variation  with  ice  displacement  were  unsuccessful,,  because  of  in¬ 
sufficient  resolution  of  the  ice  drift  track  by  satellite  navigation 
fixes.  More  quantitative  comparison  of  the  lower  frequency  transient 
currents  will  require  current  measurements  at  more  depths,  with  closer 
vertical  spacing  between  current  meters. 
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rents  first  advanced  by  Ekman,  and  »'ith  a  recent  theory  by  Pollard. 

The  inertial  currents  have  constant  phase  and  amplitude  which  only 
change  during  intervals  of  ice  displacement,  suggesting  the  local  gener¬ 
ation  of  these  currents  in  the  upper  levels  of  the  ocean. 

The  kinetic  energy  power  spectrum  derived  from  the  current  meter 
time  series  shows  an  even  fall-off  of  energy  proportional  to  the  -5/3 
power  of  frequency.  This  fall-off  is  interrupted  by  an  energy  peak  at 
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